Coexistence of bulk superconductivity and charge density wave in Cu z ZrTe3 



O 

(N 

G 



o 
o 

I 

5h 
Oh! 

c3 



c 

o 
o 



> 

o 
in 

O 



X 



Xiangdc ZhuQ Hechang Lei, and C. Pctrovic 
Condensed Matter Physics and Materials Science Department, 
Brookhaven National Laboratory, Upton, New York 11973, USA 
(Dated: June 27, 2011) 

We report coexistence of bulk superconductivity and charge density wave (CDW) with super- 
conducting critical temperature T c = 3.8 K in Cu intercalated quasi-two-dimensional crystals of 
ZrTe3. The Cu intercalation results in the expansion of the unit cell orthogonal to the Zr-Zr metal 
chains (b - axis) and partial filling of CDW energy bandgap without obvious shift of CDW transition 
temperature, b - axis resistivity pt is not related to CDW, and its dominant scattering mechanism 
for both ZrTe3 and Cuo.osZrTe3 is the electron - electron Umklapp scattering. 

PACS numbers: 71.45.Lr, 74.25. F-, 74.70.-b 



Peierls instability results in the periodic modulation 
of electron density in solids and generates an energy 
bandgap at the Fermi surface (FS), which leads to in- 
sulating behavior in one-dimensional (ID) mctal[l[. It 
can be observed as a charge density wave (CDW) su- 
perstructure in many low dimensional materials @- Su- 
perconductivity (SC), another instability of FS, will also 
generate a bandgap at FS, which results in a zero resis- 
tance. There is a view that SC and CDW orders can 
only compete with each other. NbSe3, whose structure 
is composed of infinite NbSe6 trigonal prismatic chains 
along the monoclinic b - axis, indures two CDW transi- 
tion with CDW transition temperature (Tcdw) = 145 
K (Ti) and 59 K (T2). For example, High pressure con- 
tinuously suppresses the second CDW transition and in- 
duces SC in quasi ID NbSe3 0; The suppression of charge 
order in Bai_ x A' x Bi03 via doping brings out SC with 
highest T c ~ 30 K0; Recently, both Cu intercalation^ 
and high pressure [6j in TiSe2 continuously suppress the 
CDW order, and induce SC with a dome-like phase dia- 
gram, which strongly supports the competing view. How- 
ever, CDW and SC can coexist with each other, like in 
2H-NbSe 2 with Tcdw = 33 K and T c = 7.2 K@. Angle- 
resolved photoemission spectroscopy (ARPES) results of 
2H-NbSe2 show that maximized SC at points in momen- 
tum space are directly connected by the CDW ordering 
vector, demonstrating that charge order can boost SC 
in electron-phonon coupled system Q, which is in direct 
contrast to the previous view. 

The interplay between CDW and SC is reminiscent of 
SC in heavy fermion and the cuprate oxide phase diagram 
when the magnetic order is tuned by pressure (P) or dop- 
ing to T -> Oji^. Superconducting mechanism in such 
electronic systems is likely to be mediated by the mag- 



netic fluctuations [111. Il2j. Similarly, the dome-like struc- 



ture of T c (x) and the pairing mechanism in Cu x TiSc2 is 
argued to stem from the type of quantum criticality re- 
lated to fluctuations in CDW order [lj|. Quantum phase 
transition is proposed to drive CDW order into a quan- 
tum nematic phase (13. Il5j|. On the other hand, the pos- 
itive slope of the dome could be explained via the shift 



of the Fermi level caused by Cu intercalation and the 
negative slope by the enhanced scattering for the high 
Cu concentration [3]. This is also supported by the evi- 
dence for a single gap s-wave order parameter, implying 
conventional mechanism despite the dome-like evolution 
of T c (x)[l3|- However, a dome-like structure of T C (P) 
was also discovered using P as a tuning parameter and 
suggesting that impurity effects may not be responsible 
for the dome closure @. 

Chalcogcnidc superconductors represent a weak cou- 
pling side of the smectic order, akin to stripe order in 
cuprates[I3|. There is a mounting evidence that in such 
systems CDW states transform into Fermi liquid through 
an intermediate phase [l8j|. Therefore it is of interest to 
study the melting of CDW order parameter and possible 
nematic phases in a superconductor with highly related 
and tunable two dimensional electronic system. 

ZrTe3, as an trichalcogenide (MX3, where M is the 
IV- VI transition-metal, and X is S, Se, Te), accommo- 
dates the quasi 2D ZrSe3 type structure (shown in Fig. 
1 (a))[H HtJ, rather than the quasi ID NbSe 3 type 
structure. The crystallographic b - axis defines the di- 
rection of two sets of ZrTcg trigonal prismatic chains 
with the shortest Zr-Zr distance. ZrTe3 endures a CDW 
transition at Tcdw ~ 63 K with a CDW wave vector 

1 = (&,o,j) 

anomaly on a - 



21l . |22j, and show corresponding hump 
axis resistivity (p a ) and c - axis resistiv- 



ity (p c ) around Tcrw[2l|. The FS of ZrTe3 contains one 
quasi ID electron like sheet (pocket) of states of 5p Te 
chain origin that is responsible for CDW distortion and 
pseudogap formation at T > TcdivUS SQ- ZrTe 3 
exhibits filamentary SC below 2 K with no diamagnctic 
response in H = 10 Oe(2(§. The CDW and SC in ZrTe 3 
exhibit rich and so called "competing" interplay under 
high pressure. With increasing P, the Tcdw initially 
increases below 2 GPa and then decreases up to the van- 
ishing point T C dw ~ 40 K at ~5 GPa[27j|; the T c of 
filamentary SC decreases and vanish at P ~ 0.5 GPa, 
then a reentrant SC emerges at ~5 GPa with T c rising 
monotonously up to 11 GPa (T c - 4.5 K at 11 GPa). 
Keeping in mind that the FS of NbSe3 mainly origi- 
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FIG. 1. (a)The crystal structure of ZrTe3 with blue (brown) 
symbols for Te (Zr) atoms. Comparison of the observed 
(crosses) and calculated (solid line) powder X ray patterns 
of ZrTe3 (b) and Cuo.osZrTe3 (c). The bars (black) corre- 
spond to the Bragg reflections and the lowest solid line is the 
difference between the observed and the calculated patterns. 



nates from the Nb 3d electrons and dome-like structure 
of T C (P) in TiSe2, ZrTe3 provides a different opportunity 
to study the interplay between CDW and SC. Although 
intercalation attempts have been done in alkali atoms in- 
tercalated into ZrSe^[2Sl and Cu intercalated ZrTe3 via 
electrochemical method [291] . no SC was discovered. 

In this letter, we report coexistence of bulk SC and 
CDW in Cu intercalated ZrTe3. The Cu intercalation 
results in the expansion of the unit cell orthogonal to 6 - 
axis and partial filling of CDW bandgap without obvious 
shift of Tc dw ■ Anisotropic parameters of the supercon- 
ducting state are presented, b - axis resistivity (p ) in the 
normal state is not related to CDW, and its dominant 
scattering mechanism for both ZrTc3 and Cuo.osZrTes is 
the electron - electron Umklapp scattering. 

Single crystals of Cu .o5ZrTe3 and ZrTe3 with typical 
size of 1x3x0.06 mm 3 (inset of Fig. 1(c)) elongated along 
the b - axis were grown via iodine vapor transport method 
from pure elements sealed in an evacuated quartz tube. 
After 2 days pre-reacting at 973 K, the furnace gradient 
was kept between 1023 K and 923 K. Elemental analy- 
sis was performed by using energy-dispersive X-ray spec- 
troscopy in an JEOL JSM-6500 scanning electron micro- 
scope. The Cu content is determined as the average of the 
different points on several crystals. Resistivity and mag- 
netization were measured in Quantum Design PPMS-9 
and MPMS-XL-5, respectively. The observed (Cu Ka 
radiation of Rigaku Miniflex) and calculated [3(J (Rietica 
software) powder X ray diffraction (XRD) patterns for 
ZrTe3 and Cuo.osZrTes are shown in Fig. 1 (b) and (c), 
respectively. Both materials can be indexed to the ZrTe3, 
indicating that the Cu intercalation does not change the 
crystal structure. The unit cell refinement yields lat- 
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FIG. 2. (a) p b (T) of Cu o .05ZrTe 3 near T c . Inset shows 
the M-H loop measured at T = 2 K for H_Lc. (b) (C/T) 
of Cuo.osZrTes as a function of T 2 for H = (solid circles) 
and H||c = 50 kOe (open circles). Inset shows temperature 
dependence of x'ac- See text for details. 



tice parameters a ~ 0.586(3) nm, b = 0.392(7) nm, c 
= 1.009(5) nm and /3= 97.75(1)° for ZrTe 3 , whereas for 
Cu .o5ZrTe 3 we obtained a = 0.588(2) nm, b = 0.392(8) 
nm, c = 1.011(1) nm, and j3= 97.75(1)°. The Cu inter- 
calation leads to slight expansion of the a and c - axis 
parameters. It should be noted that Zr - Zr distances 
along the b - axis chains are unchanged [20j. 

The onset ofSCatT = 3.8K can be observed in pi,, 
and the magnetization measurements confirm the SC in 
Cuo.osZrTes (Fig. 2(a)). The shape of the magnetic hys- 
teresis loop of is typical type - 2 inhomogencous super- 
conductors with some electromagnetic granularity (inset 
of Fig. 2(a)). This is in agreement with the transition 
from the temperature dependence of the real part of ac 
magnetic susceptibility (x'ac) °f Cuo.osZrTes (inset of Fig. 
2(b)), and is similar to data in some high T c cuprates 
and SmFeAs0 .85Fo.i 5 [3l], III. Although the x'ac does 
not saturate down to 1.8 K, the large estimated super- 
conducting volume fraction at T = 1.8 K (4irx' ac ~ 94%) 
reveals the bulk nature of SC in Cuo.osZrTes, as opposed 
to the pure material [26j|. 

In order to obain the electron-phonon coupling con- 
stant A , specific heat (C) of Cuo.osZrTes was measured. 
Figure 2(b) depicts the specific heat divided by T (C/T) 
as a function of T 2 of Cuo.osZrTe3 in H = kOe and 
H||c = 50 kOe. Superconducting transition is suppressed 
in 50 kOe and from its linear fit (solid line depicted in 
Fig. 2 (b)) C/T = 7 + /3T 2 ( 7 T is electron contribution, 
and /3T 3 is phonon contribution), we obtain 7= 2.64 ± 
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FIG. 3. Field dependence of magnetization of Cuo.osZrTe3 for 
H_Lc (a) and H||c (b). (c) Temperature dependence of H c i for 
H_Lc and H||c. The dashed lines are the fitted lines. Pb(T) of 
Cuo.osZrTes for H±c (d) and H||c (e). (f) H c2 (T) onset and 
H C 2(T) midpoint (see text) of Cuo.osZrTe3. 



0.12 mJ/ molK 2 and f3 = 1.21 ± 0.01 mJ/molK 4 . A Dc- 
bye temperature 6,d = 186(1) K can be estimated using 
Qd = [ i.944xU) 6 xjV ]i/3 where N ig thc number of atoms 

per formula unit. From the McMillan formula j33j 



A = 



1.04 



1.04 + ln(^f^)(l - 0.62/i*), 



(1) 



we estimate the A ~0.68 by assuming fi*= 0.13 for the 
Coulomb pseudopotential. This is a typical value of 
an intermediate coupling BCS superconductor. A spe- 
cific heat jump (AC/7T C ) can be observed around T 
= 4 K (marked by an arrow) confirming bulk SC. The 
dashed curve depicted in Fig. 2(b) is the calculated re- 
sult of the isotropic BCS gap with 2A/k^T c = 3.53 and 
AC/7T C =1.43. Although the C data in the supercon- 
ducting state is not good due to the inhomogeneity, a 
AC/7T C larger than 1.43 can be estimated, suggesting 
the intermediate or strong e — p coupling. 

Thc anisotropic properties of Cuo.osZrTea in the super- 
conducting state are investigated. Figures 3(a) and (b) 
show field dependence of magnetization of Cuo.osZrTc3 
measured at different temperatures for H_Lc and H||c. 
For H||c, the H cl is obtained using the demagnetization 
correction [34|. The H c i(T) shows almost linear realtion 
rather than H cl = H c i(0)[l-(T/T c ) 2 ] temperature depen- 
dence for both field orientation (Fig. 3(c), dashed lines). 
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FIG. 4. Temperature dependence of p a and p^ for ZrTe3 and 
Cuo.osZrTes. (a) Temperature dependence of p a (i) for ZrTe3 
and Cu .o5ZrTe 3 ; the peak (T = 57 K) and the valley (T = 
68 K) of CDW anomaly in ZrTe3 are marked with two dashed 
green lines, (b) Logarithmic plots of temperature dependence 
for Pb(i) with fitting results (dashed green lines) for the Umk- 
lapp scattering process (see text). 



Thus, we estimate H cl (0) = -0.693^-T| T=Tc from 
the Werthamer-Helfand-Hohenberg (WHH) equation[35| 
and we obtain 8.3(1) Oe (H_Lc) and 27.7(4) Oe (H||c). 
Figures 3(d) and (e) show p&(T) of Cuo.osZrTes mea- 
sured in different fields for H_Lc and H||c. Figure 3(f) 
plots the temperature dependence of H°2 set (onset) and 
H™ d (midpoint, 50% of p in the normal state). H C 2(0) 
are estimated from thc linear H C 2-T relation with WHH 
equation. Due to small variance in the in plane resis- 
tivity values plj. we assume that coherence lengths arc 
isotropic in the a - b plane (£ Q ~ £ fc = £ a6 ) and we 
estimate^ coherence leghts using H° 2 (0) = $ /(27r£ 2 fc ) 
and i?c2 b (0) = ^o/(27r^a6^c) where $ is the flux quan- 
tum. The Ginzburg - Landau (GL) parameters Kj and 
penetration depths are estimated from H l c2 (0)/H l cl (0) = 
2K 2 /lnftj (i = ab, c) and K ab (0) = A ao (0)/£ c (0), respec- 
tively. Finally, from H c (0) = H^(0)/y/2n ab (0), we es- 
timate thc thermodynamic critical field. The anisotropy 
in obtained superconducting parameters can be approx- 
imately estimated from the anisotropic GL relation: 
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where T = -^r- = — with the effective electron mass 

m ab Po-b 

m* along the i direction. The obtained mass tensor 
anisotropy r ~ 7 is comparable with the experimental 
resistivity ratio ~ 10 of ZrTe3. [2l| This points to small 
anisotropy in the gap function A(k^).[23| We list the ob- 
tained superconducting parameters in Table 1. 

The in-plane anisotropic electronic transport proper- 
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TABLE I. Superconducting State Parameters for Cuo.osZrTes. 



Cuo.o5ZrTe 3 


flj,(0) (kOe) Fji(O) (Oe) 


H c (0) (Oe) 


«i(0) 


&(0) (nm 


) Ai(0) (nm) 


, onset 
i = ab 

mid 


32(4) 
27(3) 


8.3(1) 


240(30) 
225(27) 


93(14) 
85(11) 


16(1.2) 
18(1.8) 


605(210) 
570(184) 


onset 

l = c 

mid 


13(2) 
10(2) 


27.7(4) 


330(50) 
295(50) 


28(5) 
24(5) 


6.5(1.3) 
6.7(1.3) 


1490(420) 
1539(477) 



ties in the normal state of Cuo.05ZrTe3 and ZrTe3 are 
studied. Fig. 4 shows the temperature dependence of 
p a and pb for ZrTe3 and Cu^ZrTea, which indicates the 
absolute value of resistivity is quasi 2D for both of them. 
Since the residual resistivity (po) is related to the tem- 
perature independent defect scattering contribution, re- 
sistivity subtracted by po (p(i) = p - po) can provide infor- 
mation about CDW or scattering mechanism. Compared 
with ZrTe3 (Fig. 4 (inset, a)), a weaker anomaly below 
68 K in p a (i) for Cu .o5ZrTc3 can be observed without a 
obvious shift. This is in contrast to shift of anomaly in 
ZrTe3 under high pressurej^ and lT-TiSe2 by Cu inter- 
calation or under high pressure [H, El]. Since the resistiv- 
ity anomaly (shown in Fig. 4 (inset, (a)) comes from of 
the reduction of DOS at FS due to the CDW bandgap, 
the weaker anomaly reflects that the CDW bandgap is 
partially filled via Cu intercalation. Meanwhile, the Cu 
intercalation has no significant effect on p/,(i) (shown in 
Fig. 4 (inset, b)), even though the absolute values of re- 
sistivity arc quasi 2D[2l[. The scattering along the b - 
axis in the intercalated and the pure material are not re- 
lated to the CDW. Linear relation of Pb(i) shows that the 
dominant scattering channel along Zr-Zr metal chains in 
both materials is electron-electron Umklapp scattering 
(Fig. 4 (inset, b)), as expected for quasi one dimen- 
sional chains [381] . The Umklapp process will result in a 
p(T) = po + AT™ (where 2 < n < 3) for k B T < (0.1 
- 0.3) \S\, and a linear relation at high temperature for 
k^T > 1 5 1 (3^1 . The parameter A is a constant and the \S\ 
is the typical inter-chain interaction energy. The power 
law temperature dependence is observed for T < 25 K 
with n = 2.98 ± 0.03 and n = 2.70 ± 0.03 for ZrTe 3 and 
Cuo.osZrTes, respectively (Fig. 4 (inset, b)). The esti- 
mated \S\ range is (100 — 250) K, which is consistent 
with the observed linear p&(T) above 180 K (Fig. 4). 

Generally speaking, anisotropic lattice expansion via 
Cu intercalation is different from the isotropic high pres- 
sure in ZrTc3. It reminds us that uniaxial pressure in 
NbSe3 results in different behavior with isotropic high 
pressures (3^, and that most CDW bearing materials are 
anisotropic low dimensional system. Thus, the isotropic 
high pressure provides a combined effect in these ma- 
terials, rather than a simple way of band widening as 
expected. It appears that ZrTe3 is an good candidate for 
the systematic study of the interplay and of CDW and 
SC, and further spectroscopic measurements would be of 



interest in order to unfold the details of the phonon mode 
softening and CDW melting as a function of Cu intercala- 
tion. Since SC is induced in another layered tritclluridcs 
TbTe3 under high pressure [Z(|, it is worthwhile to search 
SC by doping or intercalation in the family of CDW bear- 
ing layered rare earth tritellurides as well. 

In summary, we have shown that Cu intercalation in 
ZrTe3 results in the expansion of the lattice parameters 
in the a - c direction (orthogonal to Zr - Zr metal chains) 
and partial filling of the CDW energy gap. Bulk SC 
with T c = 3.8 K is discovered in Cuo.osZrTes, and its 
superconducting parameters are given. Along the Zr-Zr 
metal chains, electrical transport is not related to CDW, 
and the dominant dominant scattering mechanism is the 
electron - electron Umklapp scattering. 
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